INVERSE PROBLEMS IN TV-BODY SCATTERING 
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Scattering theory is the analysis of the motion of several interacting particles. 
Inverse problems in scattering theory seek the answer to the question: can one 
determine the interactions between particles by a scattering experiment? That is, 
if one shoots a number of particles at each other in an accelarator and observes the 
outcome, can one find out how the particles interact? In this note we attempt to 
explain this problem and the existing results, and we also sketch the proof of the 
extension of our recent three-body result to the many-body case: this requires only 
minor modifications. We also indicate how to extend the inverse result of [SO] to 
the many-body case. 

Before stating many-body results, we recall the simpler two-body setting, where 
(after removal of the center of mass) one studies the Hamiltonian A + V on M" for 
some n, where A is the positive Laplacian on R", and V is a function on E" that 
decays at infinity. Thus, V can be considered a perturbation of A, and hence a 
relatively simple analysis of scattering is possible. We recall that at high energies 
V can be considered a perturbation in an even stronger sense. Roughly, at high 
energies, V is not only relatively compact, but also relatively small (compared to A), 
so one can use the Born approximation to study inverse scattering. In particular, it 
is easy to recover V from the high energy (A — > -l-oo) asymptotics of the scattering 
matrix S{X) - this is the object that captures the relationship between the initial 
{t — > — oo) and final {t — > -|-oo) state of the particles of energy A, and which is 
described below. If one only knows the scattering matrix at a fixed energy A, it is 
still possible to recover V, at least for compactly supported, or indeed exponentially 
decaying, potentials EOl IH ■ For short -range potentials that are symbols, it 
is possible to recover the asymptotics of the potential at infinity using fixed energy 
information only, see For Schwartz potentials V (whose asymptotics is thus 

trivial) there are counterexamples to the fixed energy problem [H], namely there 
are 'transparent' potentials, but it would be interesting to know whether the S- 
matrix given in a finite interval determines the potential, since this would have 
direct applications in many-body scattering, see Theorem ^ 

In fact, there are many inverse problems that one can study. If there are more 
than two particles, they may have a number of initial (or incoming) and final (or 
outgoing) configurations. For example, one can shoot N electrons at each other, 
let them interact, and observe the outcome. Since the electrons repel each other, 
we will have N (asymptotically) free particles as the outcome. This is A^-cluster to 
A^-cluster scattering: there were N clusters of particles (each consisting of a single 
electron) both initially and finally. 
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In a physically much more relevant way, one may shoot an electron at a target 
which is a composite particle, for example an atom or ion. If the energy of the 
electron is high enough, the ion may break up e.g. a helium ion may break up 
into the nucleus and an electron. If the atom/ion breaks up into k parts, this is 
two-cluster to fc + 1-cluster scattering: initially there were two clusters, namely 
the electron and the ion, and finally there are fc + 1: the fc parts of the atom/ion 
and the initial electron. For the helium ion, if the ion breaks up into the nucleus 
and an electron, fc = 2, so this is two-cluster to three-cluster scattering. Or, the 
atom/ion may not break up, though it may become excited, and eventually one 
will have an atom/ion and an electron, just as initially. This is two-cluster to two- 
cluster scattering. If in addition the ion is in the same final and initial states, the 
problem may be instead considered as a two-body problem, with the ion and the 
electron comprising the two (now composite) particles; indeed, this is the usual 
way one thinks of such an experiment physically, at least at low energies. A natural 
question is then whether this is an accurate description, i.e. whether the more 
complete three-particle analysis gives the same predictions as the simplified two- 
particle approach. 

Another classification of inverse problems is by the energy range in which the 
scattering data are known. High energy problems have been extensively studied 
by Enss and Weder Novikov [21] and Wang [201 EHI- Namely, as mentioned 

above, at high energies the potential (the interaction) is small compared to the 
Laplacian (the kinetic energy), in an appropriate sense. Thus, one can use the 
Born approximation to study inverse problems and recover information about the 
interactions. In this note we discuss fixed energy problems, and problems where a 
component of the scattering matrix is known in a bounded interval of energies. 

We first recall that the actual iV-body Hamiltonian for d-dimensional particles 
takes the following form: 



where qi € K is the position and is the mass of particle i, while Vij is the 
interaction between particles i and j, and h is Planck's constant. Note that H is an 
operator on (functions on) M.^'^. The quantum iV-body problem is then the analysis 
of the Schrodingcr equation ihdtu — Hu, i.e. of the propagator e~'^*, as t — > ztoo, or 
equivalcntly, by taking the Fourier transform, of tempered distributional solutions 
of {H - \)u = 0. 

It is convenient to generalize the framework, introducing some notation due to 
Agmon and Sigal; see P| for a detailed discussion of the setup. We work on a 
vector space Xq endowed with a translation invariant metric g; by introduction 
of an orthonormal basis we may of course identify it with M", as we usually do 
below. Let A > be the positive Laplacian of g. We are also given a finite 
collection X = {Xa : a e /} of linear subspaces Xa of Xq — W\ called the collision 
planes. It is convenient to assume that X is closed under intersections, and includes 
Xq = M" and X^ = {0}. We let X" = X^ be the orthocomplement of Xa in M", 
so M" = Xa X"". We write the corresponding coordinates as {xa, x"'), and denote 
the projection to X"^ by tt". A many-body Hamiltonian is an operator of the form 



(1) 




= A + ^(^°)*K 



a 
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where Va is a real valued function on X" in a certain class, for example Va is a 
symbol on X"- of negative order: Va E S^''{X°-), p > 0. We drop the pull-back 
notation from now on and write iJ = A + Va ■ Since here we are interested in 
spectral theory, and Vq is a function on a point, hence a constant, it is convenient 
to assume that Vq = 0: otherwise it would simply shift the spectral parameter. 

In the standard example, Q, the coUison planes are Xij = {qi = qj}, as well as 
their intersections. Thus, at Xij particles i and j are at the same place, hence the 
name 'collision plane'. 

The main feature of many-body problems is that even if Va decays at infinity on 
X", it does not decay at infinity in K." since it is a constant along Xa as well as 
along its translates Xa + {x""}, a;" G X'^, so it does not decay if we go to infinity, 
say, along Xa- 

In the two-body problem one actually has H = A^j^ ^j + Vuiqi — (72), i-e. V12 still 
does not decay at infinity, e.g. if one keeps qi — 172 but lets gi — > cx). However, one 
can easily remove the center of mass by performing a Fourier transform along X12, 
and then one reduces the problem to the study of a Hamiltonian on X^^ of the 
form A + V12 where V12 decays at infinity (we are working on X^"^!). This reduced 
Hamiltonian can be considered as a perturbation of A, hence its analysis is rather 
simple. 

The center of mass can also be removed in any actual many-body problem, but 
one still obtains a Hamiltonian with non-decaying potentials as before. 

We also need to introduce subsystem Hamiltonians to describe bound states such 
as atoms and ions. We order the clusters by inclusion of the X"", so 

X" CX'' iSa< b. 

For the collision planes the relation is thus reversed, so Xa C Xi, if and only if a > 6. 
We say that a is a fc-cluster if the maximum number of elements in a chain with 
minimal element a is k. Thus, if Xa is such that Xh C Xa implies b = a or b = 1, 
then a is a 2-cluster, while 1 is a 1-cluster. If is an A''-cluster, we call H an TV-body 
Hamiltonian. It is easy to see that this agrees with the usual terminology for the 
physical iV-body Hamiltonian with center of mass removed. For example, in a 
three-body problem one gets that XaCi Xb = {0} if a 7^ 5, a, 6 ^ 0, 1. 

Corresponding to each cluster a we introduce the cluster Hamiltonian i?" as an 
operator on L'^{X'^) given by 

(2) H'' = Ax^+Y.Vb = Ax^+V'', 

b<a 

Ax" being the Laplacian of the induced metric on X". Thus, if is a A''-body 
Hamiltonian and a is a fc-cluster, then H'^ is a (A^ -I- 1 — fc)-body Hamiltonian. The 
eigenf unctions of H°- (also called bound states) play an important role in many- 
body scattering. We denote eigenvalues and normalized eigenfunctions of i?" 
by and ^'a, so (i/° — ea)'4'a = 0. The pair (a, ipa) is called a channel. We remark 
that by a result of Froese and Herbst, [H], Ca < for all a (there are no positive 
eigenvalues). Moreover, speCpp{H°') is bounded below since i?" differs from A by 
a bounded operator. Note that X^ = {0}, H'^ = 0, so the unique eigenvalue of 7?° 
is 0. 
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The eigenvalues of H'' can be used to define the set of thresholds of H"^. Namely, 
we let 

(3) Aa^Ub<aSpCCpp{H'') 

be the set of thresholds of ° , and we also let 

(4) a;, = A, U speCppiH") = Ub<a speCpp(i7''). 

Thus, e Aq for a 7^ and A^ C (— cx),0]. It follows from the Mourre theory (see 
e.g. [3123) that Aq is closed, countable, and speCpp{H°-) can only accumulate at 
Aa- In addition, eigenfunctions of i?" with eigenvalues which are not thresholds 
are necessarily Schwartz functions on (in fact, they decay exponentially, see 
0). We write 

for the intercluster interaction. Thus la is a function on M". 
The point of intruducing iJ" and Ha = Ax„ + H" is that near 

A^a,rcg A^a \ A^ajSing; 

where ATa.smg = Ub>aXb, 

Ha is a good model for H since la = H — Ha decays at infinity along Xa.rcg- (Here 
we are using that if 6 ^ a, then the collision plane Xc = Xi, D Xa satisfies c > a, 
i.e. X, C Xa.) 

If we consider the two-cluster to two-cluster problem as a two-body problem by 
regarding the atom/ion as a single particle, we essentially work on Xa by projecting 
to a bound state 4'a of and using the effective interaction 

14.cff= / /alV-aprfa;", 

which is a function on Xa- For various reasons, including the results in high energy 
settings, it is clear that we can only hope to recover Va,cS (rather than la) by study- 
ing the corresponding (i.e. state a to itself) two-cluster to two-cluster scattering 
matrix. 

We can now describe the S-matrices following [2] ■ We warn the reader that the 
normalization this gives is geometric as in jl5| . and is different from the standard 
wave-operator normalization. The difference is essentially given by pull-back by 
an antipodal map on the sphere. Thus, free particles propagate in straight lines, 
i.e. they leave in the opposite direction from where they came, so in the geometric 
formulation, in the absence of interactions, the outgoing data are supported in the 
image of the support of the incoming data under the antipodal map. An incoming 
wave of energy A in the channel a is a function of the form 

Ua - = e~'^'^~'^°l^"l \xa\ 5q _(-^)V'a(a;'^) + fastcr decaying terms, 

■ \Xa\ 

and an outgoing wave has the form 

Ua + = e'^-^^^°l^"l \xa\ ^ 5q +(T^)V'a(a;°) -I- faster decaying terms, 

i.e. the sign of the phase has changed. Here ga,± may be taken e.g. functions on 
§a, the unit sphere in Xa, or ideally, at least one of them may be taken . One 
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is then interested in generalized eigenfunctions of H, i.e. tempered distributions u 
on M" that solves {H — X)u = 0. A typical example is of the form 

(5) u = -{H-{X + iO))-\{H - A)m„,_); 

here the faster decaying terms may even be dropped without affecting u and ga,~ 
can be specified to be any smooth function on Sq . In general, even if the incoming 
data are in a single channel a, as in the corresponding generalized eigenfunction 
u of H will have outgoing waves in all channels. The S-matrix Sap{X) picks out the 
component in channel /3 by projection in a certain sense, see [3T]- Thus, 5^/5 (A) 
maps functions on Sq, the unit sphere in Xa, to functions on §{,, by 

SapWga- = 5/3,+ 

for u as in (|3J|. For example, the free-to-free (i.e. A^-cluster to A^-cluster in TV- 
body scattering) S-matrix Soo{X) maps functions on Sqi the unit sphere in M", to 
functions on §0j niore precisely 5oo(A) : L'^{Sq) L^(So) is bounded. Below we 

write §a,rog = Sa n Xa.rcg- 

The first many-body result we state, slightly informally, is an 'old' direct re- 
sult of the second author [30] and Hassell [3], which in turn can be considered as 
the extension of a result of Melrose and Zworski jJSj on asymptotically Euclidean 
spaces (i.e. in geometric two-body scattering). Its essential ingredient is the prop- 
agation of singularities, here meaning microlocal lack of decay at infinity, much as 
for real principal type PDEs in the usual microlocal setting, shown in the work of 
Hormander (lit ,10, and Kashiwara-Kawai _13^ for resp. analytic singularities 
in the absence of boundaries, by Melrose and Sjostrand ^lEl and Taylor [27j for 
smooth boundaries and singularities, and for analytic singularities by Sjostrand 
|24j . and by Lebeau [T^ for analytic singularities at corners. The relationship is 
that two-body scattering is a simpler version of the boundariless setting , 
as shown by the Fourier transform, while 3-body scattering is analogous to the 
problems with smooth boundaries, and more than three bodies correspond to the 
corners. 

Theorem 1. Suppose that H is a three-body Hamiltonian and the Va are Schwartz 
on X"" for all a. Then iSoo(A) is a finite sum of Fourier integral operators (FIOs) 
associated to the broken geodesic relation on §o to distance tt. Its canonical relation 
corresponds to the various collision patterns. The principal symbol of the term 
corresponding to a single collision at Xa is given by, and in turn determines, the 
2-body S-matrix of H"" at energies A' £ (0, A). 

This result presumably extends to short range symbolic potentials, using the 
same methods, though it is technically more complicated to write down the argu- 
ment in that case, and it has not been done. However, it was shown in |32| and j33| 
that in iV-body scattering, N arbitrary, the wave front relation of all S-matrices is 
given by the broken bicharacteristic relation. Thus, the location of the singularities 
of the Schwartz kernel of the S-matrix is known, but their precise form, e.g. that 
they are those of an FIO, is not known except in three-body scattering. (Also, in 
some parts of the wave front relation the singularities can be described precisely in 
the A''-body setting, as we do below in Theorem 0) 

An immediate corollary, when combined with two-body results (e.g. analyticity 
of the S-matrix in A' and the Born approximation) is the following inverse result. 
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Corollary 2. In three-body scattering, if the Va decay exponentially and dim Xa > 2 
for all a then Sqo{X) for a single value of X determines all interactions. 

This result is analogous to the recovery of cracks in a material by directing sound 
waves at it and observing the singularities of the reflected waves, except the last 
step which uses two-body results to get the potentials from the two-body S-matrices 
- and apart from the fact that the crack recovery result does not seem to exist in 
the literature! Note that in 3-body scattering (and also in A^-body scattering if the 
set of thresholds is discrete) there is only a finite number of possible collision pat- 
terns. Generalized broken bicharacteristics are curves in a compressed phase space 
described in or in a more leasurely way in (3^. A subset of this family is the 
set of broken bicharacteristics, which are piecewise bicharacteristics of the metric 
function g, i.e. integral curves of its Hamilton vector field in T*Xo — {Xo)x x (Xq)^, 
that satisfy in addition that at each break point at Xt.icg x ^o, the momenta ^, ^' 
of the two segments meeting there differ by an element of {X'')* , i.e. £, — £,' vanishes 
on (or is conormal to) Xt- This means that in a collision at Xh^rog, the energy and 
the component of the momentum tangential to X^ are conserved, similarly to geo- 
metric optics. A convenient way of making the broken bicharacteristics continuous 
is to compress the phase space, by identifying points (x,^) and {x,£,') ii x G Xt^rog 
and {X^)* . Such a description is in fact necessary for generalized broken 

bicharacteristics, but the new analysis below involves only broken bicharacteris- 
tics, so we adopt the more naive piecewise continuous curve approach and work in 
Xq X Xq . (In 3-body scattering every generalized broken bicharacteristic is such a 
broken bicharacteristic, but this is not true in general.) Note that a bicharacteristic 
segment of 5 is a curve of the form 

(6) x{t) = a; -I- 2t^, C = f , t G J, J an interval, 

i.e. is a straight line segment with constant momentum. In addition, is the 
kinetic energy, so in the free region, Xg^cg, |CP — ^1 but for particles in a bound 
state of energy we actually have -I- = A. 

In fact, in the phase space is compactified, much as in jJSl, which here es- 
sentially means a quotient by the K^-action on Xq x Xq that is simply the dilation 
in the Xq component. Taking So x Xq as a transversal to this action, the bicharac- 
teristic segments become curves on Sq x Xq . Below we freely use this identification 
without further commenting on it. Then the projection of a broken bicharacteristic 
segment to §0 is a (reparameterized) geodesic, and in 3-body scattering the total 
length of these segments for a maximally extended broken bicharacteristic is tt, 
see Figure ^ (One example, corresponding to an unbroken bicharacteristic, i.e. a 
straight line, is half of a great circle on §0i connecting antipodal points.) 

Unfortunately, this was necessarily a very brief and incomplete discussion, and 
we refer the reader to [23 12] for a much more detailed description of generalized 
broken bicharacteristics. 

While the A^-cluster to TV-cluster S-matrix is more complicated in general than 
for A'^ = 3, and in particular the FIO result is not known, it was proved in 
that the wave front relation of the S-matrices is given by the generalized broken 
bicharacteristic relation. These generalized broken bicharacteristics can be quite 
complicated in the presence of bound states; not all of them are necessarily broken 
bicharacteristics. However, it turns out that for an open dense subset of the relation 



INVERSE PROBLEMS IN Af-BODY SCATTERING 



7 




Figure 1. On the left, broken geodesies in M"\ {0}, n = 2, broken 
at the collision planes Xa, Xi, and Xc- On the right, the projection 
of broken geodesies in R" \ {0}, n — 3, emanating from the north 
pole, to the unit sphere §o, better understood as the sphere at 
infinity. The Co, Cf, are the intersection of the collision planes Xa, 
Xh with So; dimXa ~ 2, dimXi, — 1. 



corresponding to a single collision in an (N — l)-cluster b there are no generalized 
broken bicharacteristics relating these points except those with a single break at b. 

Indeed, it is easy to write down this relation, Af,; this was done, for example, 
in 301 Section 6]. Recall that when discussing canonical relations, such as Af, C 
r*§o X T*So, it is convenient to use the twisted symplectic form lo — to' where w, 
resp. uj', is the standard symplectic form on the first, resp. second, copy of T*So- 
We identify points in (y, fi) G r*So with pairs (y, /x) G §o x -^o = §0 x such that 
y ■ H = 0, and write for the orthogonal projection to the orthocomplement to 
span{?/}. Then the relation is the twisted conormal bundle of 

(7) {(y, y) ■■ y'b = -Vb, 2/' e So \ §&} 

with the zero section removed (twisting means that the sign of the dual variable of 
y' is changed); the restriction y' G So \ §b ensures that Q) is a smooth manifold. 
Its twisted conormal bundle is 

Ab = {{y,n,y',fi') : y' G So \S6, yi, = -y'f,, fn, = -^'f,, 
^^' = Cy\ {^,'f^-C{y')\ C = 

\ybV 

i.e. is parameterized by (y', y^, with y' G So \ Sf,, y^ G X'' satisfying \y'^\^ = 
1 — ly^P; s-iid ^^'b ^ \ {0}. Note that C is determined by y ■ fi — — y' ■ ^' . 

The (unbroken) geodesic of length tt on §o associated to (y' , /i') is the arc given 
by (spanjy'} + {c//' : c > 0}) n So; it is thus the great circle starting at y' in the 
direction /i'. Rather than using /xj^ as a parameter in JHJ, it is natural to use z G Sb 
and l/i'l = where z is the point at which the broken geodesic hits Sb- Thus, 
/i' — cP^z, c > 0, and an alternative description of A;, is 

Ab = {{y, n, y', /i') : 3z G Sb, c> 0, y' G So \ §b, 2/b = -y'b, 
H = -cP^z, fi' = cP^z}. 
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Thus, y e §o,rcg, z e §b, y'^ G X'^ with = 1 — and |^| parameterize Ab. 

(It is straightforward to express z in terms of {y' , /i') as well, and show that the 
change of variables (y',^') {y' iz) is smooth.) The reader may find it helpful to 
study the parameterizations (|11|I - H12() below. 

We also Ao be the free bicharacteristic relation, which is thus the twisted conor- 
mal bundle oi y = —y' . 

Let 

(10) 

K = {iy,f^,y',f^') ■ ^z e §>b,rcg, O O, yf, = -yi = -cPyZ, ^' = cP^z, 

y -y', y, y' i span(z) e Xa for any a 7^ 0, h}. 
Then we have the following new result, generalizing Theorem ^ 

Theorem 3. Suppose that H is an N-body Hamiltonian and the Va are Schwartz 
on X"' for all a, dimX"^ > 2 for all a. With b an (N — l)-cluster, A^, A[, as above, 
K'f^ is an open dense subset of Ah. Moreover, 5oo(A) is a Fourier integral operator 
microlocally at A'^, and its principal symbol at (y, fi,y' , fi') € A;, is determined by, 
and in turn determines, the Schwartz kernel of the S-matrix S^{X\y^\) of (which 
is a 2-body Hamiltonian) at (j^tj, ). In particular, 5oo(A) determines S^{a) for 
cr G (0, A], hence also Vh, provided that Vt decays exponentially. 

Proof. First, openness is clear since each additional condition in the definition I|1U|) 
of A'ij (as compared to (O) gives an open subset of A^, and the intersection of these 
is open. 

To see density, for z € Sb,reg, note that Xa © span(z) has codimension dimX" — 
1 > 1 in Xq, so imposing the additional conditions z e §f,,rcg and y' ^ span{z) (B Xa 
for any a ^ 0,b certainly gives a dense subset of A;,. 

Consider y' e §o,reg \ Ua^b,o(-^a ® span(2;)). Then a point in Ab is specified by 
giving y e §0 satisfying yb = —y'f,, and The projection of the corresponding 
bicharacteristic segments to the base space, Sq: i-e. the corresponding geodesies 
connecting y' and z, resp. y and z, lie in span{y', z} n So and span{?/, z} H Sq. If 
these intersect Xa then sy' + tz £ Xa, resp. sy + tz E Xa, for some {s,t) e so 
y' & Xa ® span{z}, resp. y & Xa (B span{z}, as s ^ 0. (For s — gives z g Xa, 
contradicting z G §h,reg-) Note that y' ^ Xa (9 span{z} by assumption. It remains 
to show that the set oi y G §o,rcg with j/f, = — and y ^ Xa ® span{z} is also dense 
inside 

S" = {y g So : yb = -y't}- 
Thus, 5 is a sphere inside the affine space {—y'l,} + X^, centered at — y[,, of radius 
^ Ij/bP- The other conditions on y state that it does not lie on various afhne 
subspaces Wa, namely Xa ® spanjz}, a ^ 0,6, and Wb = Xt, and in addition 
y 7^ — y'. But the intersection of a sphere with an affine space W is a lower 
dimensional sphere, unless the whole sphere lies in W. Since y = — y' does not 
lie on any of these affine subspaces Wa, « 7^ 0, all of these intersections are lower 
dimensional spheres, so in fact the set of y satisfying y ^ Xa(Bspan{z} and y G So,rcg 
is the complement of the union of compact codimension > 1 submanifolds of S, and 
hence is dense in 5*. (A slightly different way of arguing is to say that if a point y 
lies in one of these affine spaces, then the arc between y and — y' only meets this 
space inside span{y}, for otherwise — y' would also lie in this affine space, giving the 
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density rather explicitly.) Since y ~ — y' is a single point in S*, and S has dimension 
dimX'' — 1 > 1, the density follows even if we impose y ^ —y'. 

The bicharacteristics Q of the free Hamiltonian A, at energy A, can be param- 
eterized by p = (y, /i) € r*§o by 

(11) 7p W = (a^W, m). ^{t) = + 2i?, = 2/- 

The analogous 'incoming' parameterization, 7p/(i) = (x'(t), ^'(t)), p' = {y',^'), 
simply switches the sign of y and /i, i.e. it is 

(12) ^^t) = ^ + 2tt m = -^y\ 

so as t ^ +00, y along 7^, while as i -> -00, -0^ y' along 7^,. 

The normalization of these parameterizations may seem strange; they arise by 
considering the graph of d{\/Xx ■ y), resp. d{—\/\x ■ y'), (recall that e^^'^^-v ig^ up 
to constants, the Schwartz kernel of the free Poisson operator!), expressing x and 
^ in terms of {y, fj,) and an additional parameter t, and then changing the sign of 
the fi component, corresponding to the twisting. 

Now suppose that {p,p') S A'^, p — (y, /i), p' — (y',^'). Then there is a unique 
generalized broken bicharacteristic 'connecting' p and p' , i.e. with 7 such that for 
large \t\, up to reparameterization, 7 is given by 7p and 7^,; say 7|[t.+oo) coin- 
cides with 7p and 7|(_oo,T'] coincides with 7^, , up to reparameterization. More- 
over, a generalized broken bicharacteristic (being a bicharacteristic in T*Xo,rcg) 
has a unique continuation until it hits Sc for some c ^ 0, so we may assume that 
7(T),7(T') g So, sing- Note that 70, the broken bicharacteristic with a single break 
at §b (with some parameterization), is such a bicharacteristic, so ^(t) = 7o(t — i+), 
t > T, and7(i) = 7o(<-i-), t < T. In particular, 7o(r' - i_), 7o(r- 1+) e So,sing. 
Since 70 has a unique such point, they are equal, so T' — i_ — T — t+, and 
q = 7(7") ~ Since the radial momentum, j^, is strictly increasing along 

7 away from the radial points, in particular at q, we deduce that T' = T, hence 
7 = 70 as claimed. (Recall that radial points are those where the Hamilton vector 
field is radial, which in this case means that for x e ^b.rcg, = cx for some c £ M.) 

To construct 5oo(A) microlocally, we construct the Poisson operators first. Mi- 
crolocally a good approximation for these are the Poisson operators _|_(A) of Hi,, 
which are easy to write down explicitly, see j3U| . since Hj, is 'product-type': a 2- 
body Hamiltonian H'' with the 'center of mass' kinetic energy Ax^ added back. 
For suitable microlocal cutoffs (which are thus ps.d.o.'s) 

VoA^)9 = vl+(X)g - R{X + iQ)IbVl+{X)g, 
5oo(A)y = {{H - \)T+vl_iX)rVaA^): 
see [2] or Equation (4.2)]. Now, the wave front relation of 
{{H - A)r+n,-(A))*i?(A + iO)hVl4X) 

is a subset of the generalized broken bicharacteristic relation, but it includes only 
those generalized broken bicharacteristics on which /{, is not Schwartz, i.e. which 
go through points over §a with a ^ 0,b. (See the proof of Theorem 2.8 in 
Section 10] for a similar argument.) In particular, it does not contain any points 
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in A[,, SO microlocally at A'j,, Soq{X) is given by 

(13) {{H - X)T+vlJX)rvl4X) 

But this is an FIO by exactly the same proof as in the three-body setting |3()j . 
proving the FIO statement. 

The rest of the proof is exactly as in PHI Proposition 6.4], using the explicit form 
of □ 

The other new result is in 2-cluster to 2-cluster scattering. To set it up, we first 
state a theorem of Skibsted j25| . 

Theorem. |25| The two-cluster to two-cluster S-matrices have C°° Schwartz ker- 
nels, except the conormal singularity of Saa{X) corresponding to free motion, i.e. 
at the graph of the antipodal map on . 

Thus, principal symbol calculations do not help in this inverse problem. The 
new result is the following extension of [2H| to the A^-body setting. 

Theorem 4. Suppose that a is a 2-cluster, dimXa > 2, and Vt for b < a is a 
symbol of negative order (i.e. may be long range). Suppose also that 

inf spec < inf 

(and hence inf spec is an -eigenvalue of H""), I C (eQ,inf A^) is a non-empty 
open set with sup / < inf Aq, and let 

R = 2 \J sup I — Ea. 

For any fi > dim Xa there exists 6 > such that the following holds. 

Suppose that sup |(1 + \x^\YVb{x^)\ < 5 for all b with b ^ a. Then Sa'a"{X) 
given for all X ^ I and for all bound states a' ,a" of iJ" with eai,ea" < sup/, 
determines the Fourier transform of the effective interaction Va.cS in the ball of 
radius R centered at 0. 

Remark. Here 5 depends on / only through sup /; as long as sup / is bounded away 
from inf Aq, 5 may be chosen independently of /. This dependence on sup / was 
improperly omitted from the statement of Theorem 1.1 in |28| . 

The smallness assumption on the Vb, b ^ a implies that inf spec//^ is not much 
lower than inf Aa, in particular it may be assumed to lie above sup/. 

This theorem says that if the unknown interactions are small then the effective 
interaction can be partly determined from the knowledge of all S-matrices with 
incoming and outgoing data in the cluster a in the relevant energy range. In fact, 
near-forward information suffices as in two-body scattering, where this was observed 
recently by Novikov 22 , see |2H|- Also, if one is willing to take small R and a is 
the ground state of //", it suffices to know Saa{X) to recover Va^cS- That is, one 
has 

Theorem 5. Suppose that a is a 2-cluster, dimXa > 2, and Vb for b < a is a 
symbol of negative order (i.e. may be long range). Suppose also that 

ea = inf spec //" < inf Aq. 

Let ta' be the second eigenvalue of H"" , or the bottom of its essential spectrum if 
is the only eigenvalue. Suppose also that I C {ton^a') is a non-empty open set with 
sup / < inf Aa, and let 

R = 2-\/sup/ — €„. 
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For any fi > dim Xa there exists S > such that the following holds. 

Suppose that sup |(1 + |a:;''|)^V(,(a;'')| < S for all b with b ^ a. Then Saa{^) given 
for all X I determines the Fourier transform of the effective interaction Va^cS in 
the ball of radius R centered at 0. 

In case H decay exponentially on X'^ for all b ^ a, then Va^eS decays exponen- 
tially on Xa, hence its Fourier transform is analytic, so Va^cS itself can be recovered 
from these S-matrices. 

This result should extend to higher energies, i.e. sup / < inf Aa is not expected 
to be essential. But it is hard to make R greater than 2-\/inf Aa — Eq even then. 
The reason is that our method relies on the construction of exponential solutions 
following Faddeev [5| , Calderon , Sylvester and the first author [2H and Novikov 
and Khenkin , but in the many-body setting. One thus allows complex momenta 
p G C{Xa), the complexification of Xa- We then seek solutions u of {H — X)u ~ 
of the form 

(14) u^Up^e'P-^-{^^{x'')+v), p.p^X-Co., peC{Xa), 

where v is supposed to be 'small'. Note that u° = e*'''^" V'q(2:") satisfies {A + V"' — 
A)m° = for /9 as in so {H — A)^^ = /aV'a decays at infinity since tpa decays 
in cones disjoint from Xa while la decays on Xa- (We are using here that a is a 
2-cluster for otherwise la would only decay on Xa.rcgO Writing 

p = zif + p± with \iy\ ^ 1, p± ■ ly = 0, p,v real, z G C, 

the atom/ion may break up for \p_\_\ > A^ — even if A < inf A^, i.e. where 
the conservation of energy does not allow this for real frequencies. On the other 
hand, one needs such large p±^ to recover Va^cS on larger balls. 

We briefly indicate how the break-up can happen. One considers the conjugated 
Hamiltonian 

P{p) = e-*"-^" {H - A)e""^'' ^ A + 2p ■ D^^ + p ■ p + V ~ X 

with p G C{Xa) the complex frequency. The total energy of the system in state a 
is p ■ p + (a (kinetic-|-potential energy) , and we are assuming that the total energy 
is A, so 

(15) p ■ p + ea = X, see 
so 

P(p) = A -f2p y-e„. 
A good model for P(p), under the smallness assumption on /q, is 

(16) P,(p) = A + 2p-i?,„+y"-e„. 

Taking the Fourier transform in Xa makes the invertibility of Pa (p) into a question 
on the behavior of the resolvent of H"^ = A^^a -|- y(a;"), uniformly across the 
spectrum. That is, 

Tx^Pa{p)T^l = Ax. + - (e„ - l^aP - 2p • ^a), 

acting pointwise in Its right inverse is given by 

(17) Tx^Ga{p)Txl = i?"(ea - - 2p ■ ^a), 
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with i?°((T) = {H"- — a) ^, provided we show that this makes sense. Now, we need 
to keep 

(18) Fpi^a) = ea - ICaP - 2p • e C \ [inf A„ +c^), 

for otherwise other channels also become open, i.e. the cluster may break up. In 
other words, on the range of Fp, we want to keep i?° invertible on the range of 
Id— Ca, which would reflect that only the a-bound states form open channels. We 
thus need that either Im Fp(^a) 7^ or Im i^p(^a) — but ReFp{^a) < inf A^. 
Writing p = zv + p± as above, 

ImFp(^,) = -2(Imz)i/-ea, 

so for Imz 0, i.e. for p non-real, Im i^p vanishes exactly on the hyperplane v -S^a = 
0. Moreover, 

(19) ReFpi^a) - £« + pi - ((ea)± + P±)' - ■ ^af - 2(Re z){y ■ C,), 
so if lmFp{^a) ~ 0, the condition RcFp(^a) < inf amounts to 

+ pi - ((Ca)± + P±)^ < inf Aa. 
This will not be satisfied for every unless 



(20) \pi_ \ < VinfAa-e„. 

The most unfortunate restriction of the theorem is the smallness assumption on 
the unknown interactions. While H20() ensures that the 2-cluster a may not break 
up into its subsystems, it does not rule out the existence of different channels, 
e.g. associated to other 2-clusters, at complex frequencies. The presence of such 
channels would significantly complicate the analysis. In particular, the crucial 
analyticity in z would not be clear at all. 

We proceed to sketch the proof of Theorem 0] which is completely analogous to 
the three-body result shown in [2B1. For 

(21) Ai € [|pi|2 + e„,infA,)\A;, 

let ea(Ai) be the orthogonal projection to the eigenfunctions of H"^ with eigen- 
value < Ai, and let Ea be its extension to Xq via tensoring by Idx„, so 

(22) Ea{Xi)=Y,^dx^(S>it^c.'^iw)- 

a' 

Since eigenvalues of iJ" can only accumulate at A^, Ca is finite rank. We also let 

(23) Ao =inf (A;^n (Ai,+oo)) > Ai. 

The particular choice o/Ai, provided that it is sufficiently close to inf Aa, does not 
play a major role in our arguments, so we usually simply write Ca for ea(Ai), etc. 
We restrict the region 120|l slightly further and work in the region 

(24) C{Xa)l={{z,iy,p^): Imz ^ 0, jp^p + e„ e (e„, Ao) \ A^}, 

i.e. we also assume that \p±\'^ + ea is not an eigenvalue of i/". Again, we do not 
indicate Aq explicitly in the notation. 

Since the ranges of Ea and Id —Ea play a rather different role below, for p e M 
we introduce weighted spaces that reflect this: 

Hp =(L^(A:a) ® Ranea) ® (L^(A:a) ® Ran(Id-ea)) 

c{Ll{Xa)<E>Ranea)®L^Xo), 
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with Cq considered as a bounded operator on Thus, we aUow weights on 

the range of Ea, but not on its orthocomplement. Since Ranca is finite dimensional, 
Tip is a Hilbert space with norms induced on the summands by the Lp{Xa) and 
L^(Xo) norms respectively. 

Returning to Pa [p) , we can define a right inverse Ga [p) by (|17() if we show that 
it makes sense, the only issues being the behavior of R°'{a) for real a and bounds as 
l^al ^ oo. It is straightforward to analyze these, as in |2H|, and prove the following 
proposition. 

Proposition 6. Suppose that (z,i^,p±) G C(Xa)°. The operator 

Ga{p) = - If - 2p • .):Fx^ 

is a bounded operator Tip — > 7i,. for p>0,r<0,r<p— I. It satisfies 

Paip)Gaip) = Id : Hp ^ Hp, 

Ga{p)Pa{p)^ld:np^Hp. 

It is continuous in p d C(Xa)° and analytic in z € Moreover, forp > 0, r < 

0, r < p — I, p± fixed, for any G > 0, Ga{p) is uniformly bounded in B{Ti.p,Ti.r) in 
I Imz| > G\ Rez|, and s-lim|2|^oo Ga{p) — as an operator in B{TLp,TLr), provided 
that \z\ — > oo m the region |Imz| > C|Rez|. 

Sketch of proof. [See j2Hl for the completely analogous proof in the three-body set- 
ting.] Recall that We need to analyze Fp and R°'{a). First, by the spectral mapping 
theorem, for Recr < inf speciJ"^, 

\\R°'{(^)\\b(l^,l^) < (infspeciJ"-Recr)"\ 
Now, due to the — I'faP term in l|18|) . for any fixed p there exists C > such that 

(27) ReFp{^a)<G-\UV2. 

Thus, for (j) £ C^{X*) identically 1 on a large enough ball, Ga^x^{l — 4')^Xa is 
bounded on . Hence, we only need to be concerned about what happens in a 
compact set in X*. 

For this we need two other bounds that hold by the selfadjointness of i?" and 
its spectral properties, namely 

||(Id-eJi?'^(a)||6(i2,i2) <|Rec7-Aor\ Rea<Ao, 

(28) , 
W\(j)\\b(l-^,l-) <\lnia\-\ Ima^O, 

Ao as in (|23|l . It is easy to check using H19|l that 

(29) ReFp(Ca) > + pi ^ I ImFp(ea)| > (RcF,(ea) - (e^ + 

In particular, (z, v, p_i_) G C(Xa)° implies that + p\ < inf Aq, hence 

T^lR^Fp{.)){^d~Ea)Tx^ 

is well defined. Indeed, either ReFp(^a) > + p\^ and then ImFp(^a) 7^ 0, so 
R°'[Fp{^a)) is a bounded operator on L^, or Rei^p(^a) < Ca + p\ < inf A^, so 
{ld-ea)R''{Fp{^a)) is bounded on L^. 

On the range of ea,e^ the projection to the eigenspace with eigenvalue e, R"'[a) 
is multiplication by (e — cr)~^. This is a locally integrable function of tr near e (in 
C!), so the application of R°-{Fp{.)) io u ~ J-x^u is well defined, provided that 
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at every point S^a with Fp{^a) = e, the differential of Fp : X* ^ C is surjective 
(with C considered as a 2-diniensional real manifold), i.e. i^p(^a) = e implies that 
dKe Fp(^a) and dim Fp{^a) are linearly independent. But dlmFp is nonzero, and is 
a multiple of d^v-^a), so in view of H19() . dim and dRe i^p are linearly independent 
if and only if c?(((^a)_L + 7^ 0, i.e. if and only if {^a)± 7^ ~P±- But if 

(Cq)j- = -PJ- and Fp(^a) is real, then . ^„ = 0, hence Fp(^a) = £„ + pi ^ due 
to e C(Xa)°. Now let 

Ga(p) - - If - 2p • 

both terms are well defined by the preceeding considerations when applied to func- 
tions in S{Xa; L^{X'^)). Indeed, application of Ga{p) to the range of Ea,c is the 
only issue, and there, with u = v ® ■0q'j G L'^i^a), v = Tx^v^ 

so the mapping properties of Ga (p) on Ran Ea are given by the two-body results of 
Weder jSZI Theorem 1.1], see Remark 3.4 in [2HI for its application in this context. 

For the behavior of Gq(p) as |z| ^ oo, we decompose Ga as above. On the range 
of Ea we use the analogous two-body result of Weder, while on the range of Id —Ea 
we use the uniform operator bounds on i?°((T)(Id — ea), valid uniformly away 
from [Aq, -I-oo), together with | Im_F'p(^a)| oo for almost every ^a, namely the 
such that • J/ 7^ 0, so an application of the dominated convergence theorem yields 
the desired strong convergence. □ 

Since 

Pip)Gaip)v = {ld+IaGa{p))v, V € Hp, p > 0, 

we next investigate la on Tlr, again as in [2^1 . 

Lemma 7. Suppose that > 0, and {x'')^'Vb £ L°°{X^) for all b ^ a. The 
multiplication operator la is in B{Ti.r, Tip) provided that p<r + fi, p<^, r> — /i. 

Moreover, there exists G > ( independent of Vb ) such that the norm of la as 
such an operator is bounded by C max;, sup((a;'')'^|Vf,|). 

Sketch. See |28| for details. Consider the matrix decomposition of la correspond- 
ing to the direct sum in 1)25(1 . and use the rapid decay of ipa' to analyze the terms 
involving Ea. The restrictions p < r + p < ^, resp. r > — /i arise from the 
requirements of making EalaEa, Eala{l<i—Ea), resp. (Id— Ea)IaEa bounded. □ 

Thus, perturbation theory gives 

Proposition 8. Suppose that p± satisfies \p±\'^ + €a G (e^, Ao)\A^, /i > max(p, 1), 
p>0, r<0, r<p— 1. There exists 6 > with the following property. Suppose 
that for all b ^ a, sup |(x'')^Vf,| < S. Then the operator 

(30) G(p) = Ga(p)(Id +IaGaip))-' : Hp Hr 

satisfies 

P{p)G{p) ^Id-.Tip-^ Up, 
^ ' G{p)P{p) = Id : Hp ^ Hp. 

Moreover, G{p) is a continuous function of p in C{Xa)a, and an analytic function 
of z € C \ M, and s-limi^^i^o^ G{p) = as a map Tip — *■ Tir provided that \z\ —f oo 
m \ lmz\ > ClRe^l, C > 0. 
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Since /aV'a € Hp for some p > if ^ > dimXa/2, we deduce tlie following 
corollary. 

Corollary 9. Suppose that p > dimXa/2, \p±\'^ + Cq, G (eq, Aq) \ A'^. Then 

(32) Up = u°- e'''"""G(p)4Va = e^"-"" (V-a - G(p)/a^a) 
satisfies P{p)up — 0, and 

(33) w° = e*^-^''(Id+Ga(p)/a)e-''"""Up. 

Moreover, Wp — ^ m 7i,., r < 0, r < /i — 1 — dnn^fa ^ as |z| ^ cx) m | Imzj > 
G|Rez|, G> 0. 

We now only need to connect Up to the scattering matrix. This is done by letting 
p become real. Below we consider 

(34) Imz > 0, Rez > 0. 

Proposition 10. The operator Ga{p) ■ Hp — s- Hr extends continuously (from 
C(X,)+nC(X,)° j to C(X,)+ forp> 1/2, r < -1/2, and it satisfies Pa{p)Ga{p) = 
Id, Ga{p)Pa{p) = Id, on Hp, p > 1/2. The limit Ga{v, p±, z ± iO) satisfies 

e'P-^^Gaii^,p^,z±iO)e-'P-^^ 

(35) - i?a(A + lO) + Y,^xlH{v ■ - Rep))[(|eaP - (A - e^, - ^0))-^ 

- {\ia? - (A - e,, + iQ))-^]Tx^Ea,c.', 

A = + ea- 

Sketch of proof. Write 

(36) Ga{v,P±,z) = {l<l~'Ea)Ga{v,p^,z) + EaGa{l',P±,z). 

The first term is just TxlR"'{Fp{.)){ld-Ea)Tx^, and Imi^p ^ as Imz 0. 
Moreover, in the region where the sign of Imz matters, i.e. where Rei^p > Aq, 
Im i^p > by (|29|l which is valid with Inifp in place of |ImFp| due to (|34|) . so 
the first term goes to T^IR" {\£,a\'^ + 2p • + e„ + iG)[ld-Ea)Tx^ as Imz ^ 0. 
Conjugation by e*'''^" replaces by + p, hence yielding i?a(A + iQ){ld—Ea)^ 
\ = p ■ p + e^, in the limit. 

The second term on the right hand side of (|36|l is dealt with exactly as in 
the two-body setting discussed by Weder. Namely, working on the range of i?a,e, 
^Xo.Ga{p)J^x^ is multiplication by (e — Fp{^a))^^- For each with • ^ 0, the 
limit as Imz ^ is multiplication by the limit of (e + + 2p • — e^)"^, which 
is (e + I^Qp + 2p • ± *0 — Cc)"^, with + corresponding io u ■ > ^ and — to 
• < 0. Now a simple rewriting proves the proposition; see j2Hl for details. □ 

Perturbation theory now shows that G{p) itself has a limit when p becomes real, 
provided that A = p^ + < Aq. An immediate consequence is: 

Corollary 11. Suppose that \pi_\'^ +ea G (e«,Ao)\AJ,, p > (dim ATa + l)/2. There 
exists (5 > with the following property. Suppose that for all b ^ a, sup |(x'')^Vb| < 
S. Then Up extends continuously to 

(37) {z : Imz>0}U{z>0: z^ + + e„ e (e^, Aq) \ A'}, 
with Up — u'l E Hr for all r < —1/2, and Up is analytic in z in Imz > 0. 
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The last part of the proof is completely analogous to the two-body argument, 
expressing the S-matrix as a pairing, and considering an analogous pairing for the 
complex exponential solutions. For 2-clusters a, (3, the knowledge of Saf3{X) is 
equivalent to that of a renormalization (A) with free scattering removed, see 
|28[ Section 2]. Rather than giving its detailed definition here, we simply give an 
expression for its Schwartz kernel that is derived in |28l Corollary 2.2]. For p £ Xa 
(real!), we write 

Up = {ld-R{X + iO)Ia)u% X = + e^, 

Then, for A ^ A', {x'')''Vb G > dimXa, 

(39) 



5« (A,w,w')= / W^(a;)e«''-'V^(x&)dx= / hUUx)e-''''-'^ ^J^iix'') dx; 

in particular, (A) has a continuous kernel. The analogous pairing for the com- 
plex exponential solutions is 



GaQ'(p,P+C)= Ia{x)Up{x)u^^, -^^{x) dx 

(40) 

Iaix)Up{x)^Pa.'{x'')e-''-P+'^>''- dx. 

By Corollary ini if {x'')''Vb G L°°{X'') for aU b and for some p > dimX„, then 
the integral in (|40|) converges for all p for which Up exists, and for ( e M", since 
then the real parts of the exponentials cancel, and laipaipa' G ^^(^o)i f^nd the 
same holds for Iafpa'G{p){Iaipa)- Other properties of l|4Ufl follow immediately from 
Corollary El 

Proposition 12. Suppose that \p±\'^ + G (cq, Ao) \ AJ,, p > dimXa, and Vb as 
in Corollary Then Gaa' is an analytic function of z m C \ R, and extends to 
be continuous on H37() . In addition, 

(41) lim Gaa'ip,P + = / Iatpc4we~"^''''' dx, 

provided that \z\ — > cx) m |Imz| > C|Re2:|, C > 0. 
For fixed p real and C satisfying 

A = p2 + e« = (p + C)^ + e 



a" J 



i.e. the equality of incoming and outgoing energies, we can relate Gaa" (Pi P + 0? P 
real, to the S-matrices as follows. Under our assumptions, 

(42) e'P-^-{ld+Ga{p)Ia)e-"'-'''Up ^ ^ {IA+Ra{\ + lO)Ia)Up. 

Applying (Id-|-i?a(A -f iO)/a)^^ to both sides of lH^ . and using we deduce 

that 



(43) 

-f 2^l^(Id+i?,(A + lQ)Iar^TxlH{v ■ (Ca - Kep))5\^^\2^(x~,)rX^Ea.a'IaUp 



Up = Up 
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Here we used 

(p - (A - e - iO))-i - (p - (A - e + iO))-' = -27n6p_^^^,) 

with p = l^ap; the (5 distribution here is the puU-back of the standard deha distribu- 
tion on R by the map 1-^ l^ap — (A— e). Integrating against /ae^*'^''+''''^°Va"(a;°), 
and yields 

(44) -5:-=L=/ Sl^„iX,p',p + 

i?(i^-(p'-Rep))G„„,(p,p')rfp'- 

This is an integral equation for Gaa" in terms of 5^o,//(A), ^^/^/'lA), X = + £„. 
It can be solved due to the smallness assumption on la , so we deduce the following 
result, with the same proof as in \28\ . 

Proposition 13. ( 28, Proposition 5.2]j Suppose that p^+e^ = X = (p+C)^+£q" G 
(— cxD, 0) \ A', |p_Lp + £06 (cq, Aq) \ AJj, /i > dimXa. There exists 6 > with the 
following property. 

Suppose that for all b ^ a, sup | (a;'')''Vf,| < S. Then the pairings Gaa"{p,P + C) 
are determined by the operators Sf^,^,,{X) given for all a' and a" . 

We can now prove Theorem 0] Fix a non-empty open interval / C (eQ,,inf Aq) 
with sup/ < inf Aa, and let R — 2y/supl — £„. Let ( S Xa satisfy |C| < R, and 
|:ICP + ^ A^. The last condition excludes a discrete set of values of |CP in [0, R). 
Note that i|Cp + < 0. 

We now choose Ai < inf Aq in (|21|l so that ||CP + Cq < Ai, and 

(45) /n(i|CP+eo,Ai)^0. 

Let p± ~ -C/2, so p\^ ta e (ea, Ai) \ AJ,. Let j/ e be such that \v'^ = 1, 
• C = 0. Then for p = z;/ -I- p_L, z e C, (p + C)^ - = 2/9 • ( + = 0. Having fixed 
C, p±, i^, consider the energy p^ -f when z becomes real. Since p^ = + |pj_P, 
as z varies over [0, a/Aq — p^^ — ea), the energy varies over [|pj_P -I- Eq, Aq). This 
intersects the given interval / by as Ai < Aq. 

Let zo > satisfy A = p'^ + Zq + e (eo-, Aq) \ A'. By Proposition ^1 the limit 
of the pairing G^a {p, P + C) as z ^ zq (in Im z > 0) is determined by the scattering 
matrices S^^,^,, (p^ -I- Zq + e„). Thus, there exists a non-empty open interval of these 
values zq at which Gaa{p, P + C) is determined by S^,^„ (A), A G /. 

Since the limit on any open interval in the boundary of its domain determines 
an analytic function, we deduce that knowing the S-matrix S^^,^,, in the interval / 
determines Gqq(p, C) for all ( with \(\ < R. 

Now let z — > CXD through imaginary z. By (|41|l . Gaa{p, P + C) converges to 

(46) / IaHafe-'^-^^dx= f e-'^-^^(f hli^^f dxA dXa. 

But this is the Fourier transform of /alV'aP ^2;° in Xa, evaluated at C- Hence 
Sa'a"(X), X € I, determines the Fourier transform in Xa of the 'effective interaction' 
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(47) / laH^^dx'^ 

in a ball of radius 2y'sup/ — except possibly on the spheres \\C\^ + e^^ e A^. 
However, Ia\4'a\'^ dx°' S L^(Xa), hence the Fourier transform is continuous, 
hence it is determined on the whole ball. This completes the proof of Theorems 01 
and El 

Finally, we return to the electron- ion scattering experiment, comparing its treat- 
ment as a two-body problem and as a many-body problem. If one only uses a fixed 
energy A, and the potentials are exponentially decaying, in the two-body model the 
S-matrix determines the interaction, which is Va.cS in this model (here a is the 
ground state of the ion) . No such result exists if we consider the ion as a composite 
particle, i.e. if the experiment is treated as a many-body problem, although our 
theorem shows that if the 2-chister to 2-cluster S-matrix is known in an interval of 
energies and if the interactions are weak, then we can determine Ki.cff • Indeed, in 
the two-body fixed energy inverse result one has to let ^ oo. If we consider the 
scattering experiment as a many-body problem, pj^ ^ oo causes serious complica- 
tions, namely the possibility of the break up of the ion. This seems to limit the use 
of two-body results for composite particles, such as ions. 
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